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Abstract 

We present quenched results for B meson decay constants using NRQCD b quarks and 
0(a) tadpole improved clover light quarks. For the first time, one- loop matching fac- 
tors between lattice and continuum currents are incorporated through 0(a/M) taking 
operator mixing fully into account. This includes an important 0(aa) discretization 
correction to the heavy-light axial vector curr ent. We find f B = 147(ll)(±? 2 )(9)(6)MeV 
and f B Jf B = 1.20(4)(^). 

PACS numbers: 12.38. Gc, 12.39.Hg, 13.20.He, 14.40.Ndi 



1 Introduction 

The B meson decay constant fs and the bag parameter Bb play crucial roles in analyses 
of Bq — Bq mixing phenomena and the extraction of CKM matrix elements. Experimental 
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determination of fs is severely constrained by the very small leptonic branching fraction, 
while Bb is not a direct outcome of experiments. Their calculation has proven difficult due 
to the possibility of large QCD corrections to the basic weak process. Lattice QCD provides 
a non-perturbative method for calculating the relevant hadronic matrix elements from first 
principles fl|] . This paper presents a quenched lattice calculation of fs and fs s using the non- 
relativistic QCD (NRQCD) approach to 6-quarks 0]. Much of our discussion of the systematic 
improvement of Jb, however, applies to other lattice calculations of this quantity. 

In the rest frame of the B meson, the momentum of the o-quark is of O(Aqcd)- The 
6-quark velocity is hence v « Aqcd/M m 0.1, where M is the heavy quark mass, and the non- 
relativistic approach is sensible. The heavy quark action that we employ includes relativistic 
corrections through 0((Aq CD /M) 2 ). For the light quarks we use the tree- level tadpole im- 
proved clover action which reduces 0(clAqcd) errors to 0(a (clAqcd)) 0- In the heavy-light 
axial current we investigate both 1/M and 1/M 2 corrections. The one- loop matching factors 
relating the lattice and full (relativistic) continuum QCD currents have been calculated to 
0(a/M) 0] and incorporated in this analysis. Mixing among heavy-light current operators 
and contributions from an 0(a(aAQCD)) discretization correction to the current are taken 
into account. In the nonrelativistic approach the 1/M expansion and 0(a) improvements go 
hand in hand. Previous calculations of fs using NRQCD 6-quarks have included only sub- 
sets of the above ingredients || |6], |7j, as have calculations using the relativistic formalism for 
the 6-quark |], Rl [H], 11, 12]. Here we include all of them and carry out a comprehensive 



study with improvements through 0((Aq C d/M) 2 ) and 0((iAqcd) in the action and through 
0((A QCD /M) 2 ), 0(«/(aM)), 0{a{A QCD /M)) and 0{a{aA QCD )) in the currents. 



2 Details of the Simulation 



The NRQCD heavy quark action density used in our lattice simulations is a variant of the one 
proposed in [[J, 

„ - - / a5H\ ( aH \ n + / aH \ n ( a5H\ 

where ip denotes a two-component Pauli spinor, Ho is the nonrelativistic kinetic energy oper- 
ator, 

#o = , (2) 

2M ' V ; 

and 5 H includes relativistic and finite-lattice-spacing corrections, 

5H = ?^<t-B 

2M 

"•' (V ■ E - E ■ V) - ,? r , 9 cr ■ (V x E - E x V) 
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In addition to all 1/M 2 terms we have also included the leading 1/M 3 relativistic correction as 
well as the discretization corrections appearing at the same order in the momentum expansion. 
V and A^ 2 ) are the symmetric gauge-covariant lattice derivative and laplacian, while is a 
discretized version of the continuum operator Df. The parameter n is introduced to remove 
instabilities in the heavy quark propagator due to the highest momentum modes p|. For the 



heavy quark action, we use two-leaf E fields [DJ rather than the four-leaf fields of previous 



NRQCD simulations. The tree-level coefficients are tadpole-improved by rescaling the gauge 
fields by the fourth root of the plaquette, uq = 0.87779 at (5 = 6.0 Ul4| . For the light quarks 
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we use the tree-level tadpole improved clover action (csw — u o J- 

Simulations were carried out on 102 quenched configurations on 16 3 x 48 lattices at (3 = 6.0. 
Both forward and time reversed quark propagators were calculated on each configuration in 
order to improve statistics. Error estimates are obtained using a bootstrap procedure. Meson 
correlators are calculated for 30 combinations of quark masses; five values for k and six values 
for the bare heavy quark mass. Two of the k values are selected to simulate the charm quark, 
and the remaining three are around the strange quark mass. Heavy quark masses range from 
below the bottom quark mass to about five times the bottom quark mass. Both the heavy 
and the light quark propagators were smeared at the sink and/or the source. Further details 
are given in Ref. |J. 

There are four input parameters in the simulations, the light quark mass (/fy = k u ^), the 
strange quark mass (k s ), the scale a -1 , and the heavy quark mass (M ). These need to be 
fixed via experimental input. 

To set the scale one can use a variety of quantities like M p , f n , string tension (a), or 
quarkonium S-P splittings. In quenched calculations, however, different physical quantities 
can lead to different a -1 's and one must make a choice. In principle one would like to choose an 
observable whose characteristic momenta are similar to those of the quantities being calculated 
and which has comparable quenching and discretization errors. The first choice would be to 
fix a -1 from the heavy-light spectrum itself but experimental and numerical errors will have 
to be reduced before this becomes practical. Based on the "brown muck" picture of HQET 



for heavy-light systems fLo| , however, we believe that the scales relevant for energy splittings 
in heavy-light spectroscopy and for the decay constant are determined by the light quark 
dynamics. For Jb and spectroscopy we shall therefore consider observables in the light quark 
and gluonic sector to set the scale. For clover fermions at (3 = 6.0, the values for a -1 from M p , 
M K , M N , M A , a, and / w range between 1.8 GeV and 2.0 GeV. We use a~* = 1.92(7) GeV as 
it is close to the mean, and can be extracted most reliably from our data. We also calculate 
the final values using a" 1 = 1.8 and 2.0 GeV and use the spread as the uncertainty due to the 
determination of the scale. 

We determine Ki = 0.13917(9) by linearly extrapolating the pion mass, calculated on the 
same set of configurations, to the physical value. To set k s we use the K, the K* and the 
(f), which yields the values 0.13755(13), 0.13719(25), and 0.13717(25), respectively. We quote 
central values for our B s meson decay constant using the K, and the difference between the 
K and 6 systematic error. 

The bare heavy quark mass is fixed by linearly interpolating between the meson masses with 
the three lightest M values to the experimental B meson mass. Since the NRQCD action 
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u > J -'sim, 


aM kin 


aM per t 


aM' 


a E ■ 
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aM kin 


aM pert 


aM' 


1.6 


0.422(6) 


2.17(24) 


2.08(09) 


2.11(5) 


0.483(8) 


2.21(13) 


2.14(09) 


2.17(5) 


2.0 


0.438(7) 


2.57(33) 


2.46(10) 


2.52(6) 


0.499(8) 


2.63(18) 


2.52(10) 


2.58(6) 


2.7 


0.452(8) 


3.34(58) 


3.12(11) 


3.26(8) 


0.513(8) 


3.37(29) 


3.18(11) 


3.32(8) 


4.0 


0.461(9) 


5.2(1.5) 


4.32(10) 


4.59(9) 


0.522(8) 


4.86(67) 


4.38(10) 


4.65(9) 



Table 1: Simulation energies and pseudoscalar meson masses in lattice units for light and 
strange clover quark masses. Meson masses are determined from eq.@ with A coming from 
perturbation theory (M pert ) or heavy-heavy spectroscopy (M'), and from the dispersion rela- 
tion eq.(|) (M kin ). 

omits the heavy quark rest mass, heavy-light pseudoscalar meson masses are not direct outputs 
of the simulations. Instead, correlation functions fall off with an energy E sim that is related 
to the full meson mass Mp$ and the bare 6-quark mass M as, 

M PS = A + E sim = [Z m M - E ] + E sim , (4) 

where Z m is the mass renormalization constant and E an energy shift. Z m and E Q have been 
calculated perturbatively and E sim is obtained nonperturbatively from the simulations. An 
alternate way to obtain the meson mass is to calculate energy splittings between correlators 
with and without spatial momentum, 

E sim {p) - E sim (0) = ^ + Ml m - M kin . (5) 

In Table [l] we list results for Mps of eq.(fD using perturbation theory to calculate A (denoted 
" Mpert 1 in the Table) and for M k i n of eq.@, both measured in units of a -1 . A can also be 
obtained nonperturbatively from simulations of quarkonium QQ systems, whose correlations 
fall off with an energy Ef in \ related to the full QQ meson mass as, 

Mq Q = M^ = 2A + E^. (6) 

Both MgJ and E%% can be extracted from the simulations, leading to a nonperturbative 
determination of A. The resulting masses for the heavy-light mesons are also listed in Table [TJ 
under aM'. One finds consistency among the three determinations of the meson mass, aM kin , 
aM pert and aM', in most cases, to within la. The errors for aMkin are much larger than for 
aM' since finite momentum correlators have much larger fluctuations for heavy-light mesons 
than for heavy-heavy mesons. The errors on aM pert were estimated by squaring the one-loop 
corrections and by varying q*, the scale in the coupling a. We will use the most accurately 
determined aM' to fix aM . With a _1 = 1.92(7) GeV and using linear interpolation between 
the three lightest heavy quark masses, one finds that aMo = 2.22(11) gives the correct B 
meson mass. Had one used aM^in rather than aM', the number changes to aMo = 2.28(42). 
The effect of this difference on fg is invisible. 

The results for the heavy-light spectrum, including radially excited states, P-states and 
heavy baryons, based on the same data set, are summarized in Figure 1. We find that the 
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Figure 1: Heavy- light spectrum with 6-quarks. Horizontal dashed lines bracket experimental 
values taken from the Particle Data Book. Horizontal dotted lines indicate first experimental 
observations of these states ([|T7]] (2S & \ +l / B^) and |TB[ Errors on simulation results 

include statistical plus, where appropriate, errors from extrapolation in the light quark mass 
to the physical pion. The two results for the strange B mesons correspond to fixing the strange 
quark mass from either the K (filled circles) or the <fi (open circles). For the J = 1 P states 
only one level is shown, since we are currently unable to reliably separate the physical 1 + and 
1 +/ states. 



combination, NRQCD heavy and clover light quarks, reproduces the gross features of the B 
mesons and the heavy baryons. This is a prerequisite for a reliable calculation of fs- Details 
of the spectrum calculations will be given in a separate publication [|16j. 

3 Heavy-Light Axial Currents in NRQCD and Match- 
ing to Continuum Full QCD 



The goal is to calculate the pseudoscalar decay constant f PS , defined in Euclidean space 
through the relation, 

(0\A tM \PS)= Ptl f PS . (7) 

Here is the heavy- light axial vector current in full continuum QCD. One needs to relate 
this current to the current operators jjp in lattice NRQCD. 

Heavy-light currents in full continuum QCD have the form qTh. For Aq, we use T = 7570. The 
four-component Dirac spinor for the heavy quark, h, is related to the two-component NRQCD 
heavy quark (heavy anti-quark) fields, ip (ip), via an inverse Foldy-Wouthuysen transformation: 
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h = u^ FW = u^[ Y ~ y (8) 

Through 0(1/M 2 ) this gives, 

qTh = qTQ-^qT( 1 -B)Q 

J-gr(D 2 + ^S-B-2^a-E) Q, (9) 



8M 2 



where Q = |(1 + 'Jo)^fw and a = j j, £ = diag(cr, cr). In our simulations we have included 
all tree-level contributions to the current through 0(1/M 2 ) using the terms in eq.(fj|). 

We use perturbation theory to match between the lattice currents in our simulations and those 
of continuum QCD and employ on-shell matching conditions. One starts by considering the 
process in which a heavy quark of momentum p is scattered by the heavy-light current into a 
light quark of momentum p' . Matching takes place by requiring that the amplitude for this 
process calculated in full QCD agrees with that in the effective theory through the order in a 
and 1/M that one is working in. Here we work through 0(a/M). The amplitude in full QCD 
can be expanded in terms of p/M, p'/M etc. The next step is to identify current operators 
in the effective theory that would reproduce the same terms in the amplitude. Finally, a 
one-loop mixing matrix calculation must be done within the effective theory. Details of such 
a matching calculation are given in 0. Here we summarize the main features. 



The full QCD calculation uses naive dimensional regularization (NDR) in the MS scheme. A 
gluon mass A is introduced at intermediate stages of the calculation to regulate IR divergences. 
This is permissible here since no non-Abelian vertices appear in this one-loop calculation. A 
one-loop calculation in full QCD finds, upon expanding through 0(1/M), that the following 
three current operators are required in the lattice effective theory: 

J? = -^9757o(7 "V)Q, 

^g(7-V) 7 57oQ. (10) 



7(2) 



jf^ and Jl correspond to discretized versions of the first two terms in eq.(|9|), appears 
only at one loop. After carrying out a one-loop mixing matrix calculation in lattice NRQCD, 
one ends up with the final matching relation between (Aq) in full QCD and matrix elements 
evaluated in the lattice simulations, 

(^olgc^E^J?). (11) 



Co, C\ and C2 are the three matching coefficients. Relation (|TT| ) is independent of the external 
states up to lattice artefacts. 
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An interesting outcome of the one-loop calculation is that the mixing between and 
does not vanish in the limit M — > oo. This is because of the generation of a term, at order(a a), 
which is a lattice artefact. This term can be absorbed into an 0(a (qAqcd)) discretization 



correction to J L 



(o). 

(°) . 7(°) i n r( disc ) 



Jl -> Jr+C A Jr C \ (12) 
with 

jt C) = aqh-^hsloQ- (13) 

The coefficient Ca is also calculated in ||. j| dlsc ) i s the analogue in heavy-light physics of 
the a <9 M P discretization correction to the light quark axial current (here P is the pseudoscalar 
density) of Ref. jn5 . 



The matching coefficients Cj are functions of the coupling a and of the heavy quark mass and 
one needs to specify the precise definitions of these parameters. We use cty(g*) of reference 



14 1 as our perturbative expansion parameter. The scale q* is not yet known for this matching 



calculation, however, a value for the M — ► oo limit has been calculated in reference [£0] and 
gives q* = 2.18/a. We shall extract fg for q* = 1/a and q* = n/a and use the difference 
between the two as an estimate of the a 2 uncertainty. It is also important to use a common 
definition for the heavy quark mass in the continuum and lattice calculations. This can be 
ensured in perturbation theory by employing on-shell quark mass renormalization and using 
the pole mass in both theories. After writing everything in terms of the pole mass, we convert 
to the bare lattice mass Mo, which is a well defined quantity even beyond perturbation theory 
and the mass that appears in our lattice simulations. We use one-loop perturbation theory 
for this conversion. 

Once 0(1/M 2 ) corrections to the current are included at tree level, eq.([Ll|) is modified to, 

(4>)qcd = E Cjijjp) + E (Jl), (14) 

.3=0,1,2 .7=3,4,5 

with j]p , j = 3,4,5, corresponding to the three 1/M 2 current corrections in eq.(|9]). The 
matching coefficients Cj, j = 0, 1, 2, do not include one-loop feedback from these latter 1/M 2 
current operators. Such contributions come in at 0(a/M 2 ). Eq. (Ill) is the NRQCD/Clover 
heavy-light axial current we use to extract f PS . It includes terms through 0((Aq C - £) /M) 2 ), 
0(aA QCD /M), 0(a/(aM)) and 0(a (aA QCD )). We note that powers of both l/(aM) and 
Aqcd/M are contained in our expressions. These ratios are both considered small in the 
NRQCD formalism (in practice l/(aM) is sometimes allowed to become as large as ~1). In 
particular, one does not take the limit aM — > 0. The cutoff and therefore momenta allowed 
cannot become larger than ~ M in a non-relativistic theory. 

We end this section with two brief comments. The first concerns non-perturbative effects 
that arise because we have power-law contributions to the matching coefficients that go as 
a/(aM) n . These terms are not separated out explicitly in our calculations but they are evident 
if coefficients are plotted as a function of aM. They are there to cancel unphysical ultra-violet 
pieces in the NRQCD current matrix elements, which do not match full QCD. It is possible 
for the a/(aM) n terms in, say, C to multiply non-perturbative lattice errors, oc (oAqc'd)") in 
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the matrix element ( J L ) to give contributions that look like physical (Aqc^/M)™ corrections 
to Jb- This limits our ability to extract the physical M dependence of fs- The source of these 
non-perturbative errors is a mis-match between the infra-red physics of lattice NRQCD and 
that of full QCD. After taking the improvement of eq.(|T2"D into account our is designed 

to match full QCD through 0(a a). So these errors appear at 0(a 2 (ciAqcd)) or 0((a Aqcd) 2 )- 
The leading contribution to fn from these non-perturbative lattice artefacts is then ajiaM) x 
a 2 (oAqcd) = « 3 X Aqcd I 'M or a/ \aM) x (aA QCD ) 2 = a (aA QCD ) x A QCD /M giving 
an error in the slope of 0(a 3 ) or 0(a {clAc-cd))- Both these terms are of the same order as 
higher order corrections not yet included in our matching calculations; they do not lead to any 
fundamental difficulties. Better precision can be attained in principle by going to more highly 
improved actions and current operators. Conversely we note that the situation is considerably 
worse for Wilson light quarks where discretization errors in the action mean that the currents 
cannot be improved to the level we have achieved here. Then non-perturbative lattice artefacts 
could induce an a/(aM) x (aAqcD) = ol (Aq C d/M) error in f B . 

The second comment concerns cancellation of the power-law contributions referred to 
above. For instance, our matrix elements will contain unphysical contributions that behave 
as a/(aM) and which are cancelled by the perturbative matching procedure. This might be 
a delicate operation if the terms behaving as a/(aM) in the matching coefficients and in the 
matrix elements were much larger than the physical terms of 0(Aqcd/M). In the region in 
which we work none of the perturbative coefficients of the a/(aM) terms are unduly large. No 
delicate cancellation is necessary. Furthermore the current matrix elements are clearly domi- 
nated by their physical components, in particular the matrix elements of the higher dimension 
operators, and , show no sign of significant ajiaM) 11 terms. This will be discussed 
in more detail elsewhere pi| . In our calculations power- law terms have been subtracted per- 
turbatively through 0(a/(aM)). Uncancelled contributions would come in at 0(a 2 /(aM)) 
and 0(a/(aM) 2 ). There is no indication from the behavior of our matrix elements that such 
contributions are particularly large. Any uncertainty arising from them will be covered by 
the systematic errors that we assign in the next section to higher order perturbative and 
relativistic corrections. 



4 Results 



Details of our data analysis for decay constants are described in []5], |(J and will not be repeated 
here. We start with the discussion of the size of the various terms in the lattice current defined 
in eq.(|TJ). Table ^ summarizes our results for 

/«> JmJs = — != (0| J? |PS>, (j = . . . 5), (15) 

V M PS 

evaluated at zero 3-momentum. The light quark has been extrapolated/interpolated linearly 
to Ki or k s . The matrix elements of are not shown as they are identical to at zero 3- 
momentum. The data show that around the physical B (aMo = 2.22(11)) the matrix elements 
of are ~ 12% of \ The analogous ratio for the sum of the three 1/M 2 currents (which 
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aM 




a 3 / 2 / (1) VM 


a 3 / 2 / (3) >/M 












1.6 


0.164(7) 


-0.0245(12) 


-0.0053(4) 


0.00363(15) 


-0.0050(3) 


2.0 


0.167(9) 


-0.0206(10) 


-0.0036(2) 


0.00216(10) 


-0.00325(19) 


2.7 


0.170(11) 


-0.0162(10) 


-0.00220(14) 


0.00109(6) 


-0.00181(13) 


4.0 


0.174(11) 


-0.0113(8) 


-0.00111(8) 


0.00045(3) 


-0.00082(5) 


7.0 


0.191(8) 


-0.0075(4) 


-0.00049(2) 


0.000144(8) 


-0.000287(15) 


10.0 


0.198(9) 


-0.0055(2) 


-0.000272(12) 


0.000072(5) 


-0.000145(8) 








1.6 


0.194(5) 


-0.0278(8) 


-0.00635(17) 


0.00420(11) 


-0.00576(17) 


2.0 


0.198(6) 


-0.0236(7) 


-0.00446(15) 


0.00253(7) 


-0.00375(12) 


2.7 


0.205(7) 


-0.0187(6) 


-0.00276(10) 


0.00129(4) 


-0.00210(8) 


4.0 


0.212(8) 


-0.0134(5) 


-0.00144(6) 


0.00054(2) 


-0.00098(3) 


7.0 


0.230(6) 


-0.0088(2) 


-0.000620(17) 


0.000171(5) 


-0.000341(10) 


10.0 


0.237(6) 


-0.0065(17) 


-0.000344(9) 


0.000083(3) 


-0.000171(5) 



Table 2: Tree-level matrix elements, in lattice units, is defined in eg. (|T5|) . 

partially cancel) is ~ 3%. It is therefore reasonable to assume that the neglected 1/M 3 and 
higher terms are < 1%. 

The results for a 3//2 fps a/Mp^ with the one-loop matching coefficients are summarized in Table 
3. We show results for three values of a: 0, ay(q* = l/a) and ay(q* = ir/a). Two of the 
matching coefficients, Co and C\, include a term a \n(aM)/n. In Table 3 we have used the 
dimensionless bare heavy quark mass, aM , for the argument of these logarithms. There are 
no problems with large logarithms around the physical B meson, however, if one were to work 
at large aMo, resumming the logarithms would be more appropriate. At the B meson, the one- 
loop correction decreases the decay constant by ~ 9 — 14% depending on q* of which ~ 3 — 5% 
are due to the one-loop contributions from and ■ The ~ 5% uncertainty 

coming from the variation in q* will be taken as a measure of the uncertainty coming from 
higher order perturbative corrections. This is a conservative assignment of ~ 50% of the one- 
loop correction. The perturbative matching coefficients employed in this article used the full 
1/M 2 action of eq.(0) with 8H given by eq.([|) (with, however, four-leaf E fields). Matching 
coefficients from an 0(1/M) NRQCD action (with 5H = —^m^ " ' B) are also available [|J 
and lead to a difference in fs of about ~ 1%. This comparison gives us some idea of the effect 
from 0(a/M 2 ) contributions to matching coefficients on We will assign a ~ 4% error to 
the full 0(a/M 2 ) corrections from the 1/M 2 current operators. This is the size of a/(aM ) 2 
and is more conservative than using the tree level value of the 1/M 2 corrections. Adding this 
in quadrature to the ~ 5% 0(a 2 ) and ~ 1% 0(1/M 3 ) errors gives a ~ 6% systematic error 
from higher order perturbative and relativistic corrections. The dominant discretization errors 
in the present decay constant calculation are of 0((aAg C £>) 2 ), giving an additional systematic 
error of ~ 4%. 
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ki 




aM 


tree-level 


q* = n/a 


q* = 1/a 


tree-level 


q* = n/a 


q* = 1/a 


1.6 


0.133(7) 


0.121(6) 


0.114(6) 


0.158(5) 


0.144(4) 


0.136(4) 


2.0 


0.142(8) 


0.128(7) 


0.121(7) 


0.169(5) 


0.153(5) 


0.144(5) 


2.7 


0.151(10) 


0.135(9) 


0.127(9) 


0.182(7) 


0.164(6) 


0.153(6) 


4.0 


0.162(10) 


0.144(10) 


0.135(9) 


0.196(7) 


0.176(6) 


0.164(6) 


7.0 


0.183(8) 


0.162(7) 


0.150(7) 


0.220(6) 


0.195(5) 


0.181(5) 


10.0 


0.193(9) 


0.170(8) 


0.158(8) 


0.231(6) 


0.204(6) 


0.189(5) 



Table 3: Decay matrix elements fy/M in lattice units. The first column gives tree-level 
results, the second includes renormalization constants using ay at q* = tx /a, and the third 
using ay at q* = 1/a. 







f B [MeV] 


fs. [MeV] 


a- 1 [GeV] 


aM 6 ° 


q* = n /a 


q* = 1/a 


q* = ir/a 


q* = 1/a 


1.8 

1.92 

2.0 


2.39(7) 
2.22(6) 
2.11(6) 


140(9) 
151(9) 
159(9) 


131(8) 
142(9) 
150(9) 


171(5) 
180(6) 
187(6) 


160(5) 
169(6) 
176(6) 



Table 4: Decay constants f B and f Bs in MeV. Statistical errors only are shown. aM 6 ° is the 
dimensionless bare heavy quark mass which leads to the correct B meson mass. 

Table [| shows results for fp$ at the physical 6-quark mass for several choices of the scale. For 
each scale, aM 6 ° denotes the value of the dimensionless bare heavy quark mass which gives the 
correct B meson mass. Using the scale a~* = 1.92(7) GeV and averaging over the two q*'s, 
our final estimate in the quenched approximation is then, 

f B = 147(11) (±l 2 ) (9) (6)MeV, 

f B , = 175(08)(^ )(ll)(7)(l 7 )MeV. (16) 

The first error comes from statistics plus extrapolation/interpolations in k and aM Q . Scale 
uncertainties are reflected in the second error with the upper and lower limits coming from 
a -1 = 2.0 GeV or a -1 = 1.8 GeV respectively. The third error is due to higher order relativistic 
and perturbative corrections and the fourth due to discretization corrections. The central value 
of f Bs has been calculated fixing the strange quark mass from the K meson; the fifth error 
bar comes from the difference in k s from the <j) and the K. 

Another quantity of interest is the ratio fsj ' fs- It can be determined to a greater accuracy 
than fs and f Bs separately, since the direct dependence on the lattice scale drops out, and the 
statistical errors partly cancel. We obtain this ratio for each value of the heavy mass using, 

U (fVM) K=Ks ^M~ B 

fs (fVM) K=Kl vmf: 1 } 
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aM 


tree-level 


q* = n/a 


q* = l/a 


1.6 


1.18(4)(5) 


1.18(4)(5) 


1.19(4)(4) 


2.0 


1.19(4)(4) 


1.19(4)(4) 


1.19(4)(5) 


2.7 


1.20(5)(5) 


1.20(5)(5) 


1.21(5)(5) 


4.0 


1.21(5)(5) 


1.21(5)(5) 


1.21(5)(5) 


7.0 


1.19(4)(5) 


1.20(4)(4) 


1.20(4)(4) 


10.0 


1.19(4)(4) 


1.19(4)(5) 


1.19(4)(5) 



Table 5: fsj fB, without renormalization constants (left), q* = n/a (middle), and q* = l/a. 



where M B and M Bs are the experimental meson masses. Table || shows our results. One 
sees that one-loop renormalization has negligible effect on f B Jf B . This also indicates that 
the errors from higher orders in perturbation theory are negligible. In the range studied, the 
data do not show any significant dependence on aM . Linearly interpolating to aM°, our best 
estimate of the ratio is 

^ = 1.20(4)0*). (18) 

JB 

The first error is statistical and the second error comes from the uncertainties in fixing k s . 

Eqs. ([If]) and flIB|) give the main results of this article and represent the most complete results 
to date for f B and /b„//b using NRQCD 6-quarks. It is interesting to compare them with 
other calculations that have used a relativistic formalism for the heavy 6-quark. Table [5] gives 
such a comparison with the APE §, Fermilab M, JLQCD pf, MILC |TJ and UKQCD 



12| collaborations. Details of the simulation parameters and analyses are different among 
these groups (for recent reviews see |1|). For example MILC uses the Wilson quark action, 
UKQCD and APE results are for the clover action with csw — 1, Fermilab uses tree- level 
tadpole improved clover and JLQCD uses one-loop tadpole modified clover. Only JLQCD, 
MILC and Fermilab extrapolate to a = 0, and JLQCD and Fermilab include M-dependent 
renormalization constants at O(a). They do not take into account one-loop mixing with 1/M 
current operators or the discretization correction that arises at 0(a(aAQCD))- We expect that 
at comparable (3 and M these two corrections should be of similar size in their formalism as 
we find here (~ 3 — 5%). In spite of the range of approaches and approximations used, it is, 
nevertheless, encouraging to see that very different formulations of heavy quarks give results 
for f B and f Bs that are consistent within errors. The ratio f Bs / f B , in which many of the errors 
cancel, shows a larger spread among the groups. The NRQCD number lies 1 ~ 2a higher than 
the more recent relativistic heavy quark results in Table |6| by APE, Fermilab and MILC. 

Our B meson decay constant results are still within the quenched approximation and at a 
fixed value of the lattice spacing. Calculations are underway to address these sources of un- 
certainty. The first results from nf = 2 dynamical configurations (without the 1/M 2 tree-level 
corrections) show that, within current statistical and systematic errors, unquenching effects 



are not yet visible p2| . Higher statistics and smaller dynamical quark masses are required to 
resolve this issue. Since we are using improved actions and currents we expect scaling viola- 
tions to be small in quenched calculations at (3 = 6.0. Further quenched simulations at other 
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Collaboration 


/b [MeV] 


/b. [MeV] 


/b s / /b 


This work 


147(1^) 


175(111) 


1.20(12) 


APE § 
Fermilab @ 
JLQCD (TJ 
MILC QTJ 
UKQCD |TJ 


180(32) 
164(111) 
163(16) 
153(1?°) 
160(11) 


205(35) 
185(11?) 
175(18) 
164(lf°) 

194(+g) 


1.14(8) 
1.13(11) 

1.10(6) 
1.22(4) 



Table 6: Comparison with quenched results for /b, /b s and /b s //b from recent lattice cal- 
culations using relativistic quarks. Errors include statistical and systematic errors, added in 
quadrature. 



5.7 and j3 = 6.2) have been initiated to check this. Preliminary results 

however, to make any definitive 



lattice spacings (/? 

at (3 = 5.7 on UKQCD configurations are encouraging |23 
statement requires a full analysis at all three f3 values. 



5 Summary 

This article describes a quenched calculation of B meson decay constants, fs and /b s , using 
NRQCD 6-quarks and clover light quarks. Tree-level corrections to the heavy quark action and 
the heavy-light currents have been included through 0(1/M 2 ). The full 0(a/M) matching 
between continuum QCD and lattice NRQCD currents has been incorporated together with 
an O(aa) discretization correction to the heavy- light axial current. This is the first time both 
mass-dependent matching factors and the O(aa) current correction have been included in a 



lattice determination of fs- Our final results for /b, /b s and /b s //b are given in eqs. ([Tq) 
and (|I5D . The next step in our program is to study unquenching, scaling and finite volume 
effects, in addition to improving statistics. Work is already in progress on these fronts. 
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